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ABSTRACT 

We report the discovery of a 2.1hr optical modulation in the transient source GS1826- 
24, based on two independent high time-resolution photometric observing runs. There 
is additional irregular variability on shorter timescales. The source also exhibited an 
optical burst during each observation, with peak fluxes consistent with those of the 
three X-ray bursts so far detected by BeppoSAX. We compare the low-amplitude 
variation (~ 0.06 m ) to that seen on the orbital periods of the short period X-ray 
bursters, X1636-536 and X1735-444, as well as the similarity in their non-periodic 
fluctuations. Other transient neutron star LMXBs possess short periods in the range 
3.8-7.1 hrs. However, if confirmed as the orbital, a 2.1 hr modulation would make 
GS 1826-24 unique and therefore of great interest within the context of their formation 
and evolution. 
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1 INTRODUCTION 



GS 1826-24 was discovered serendipitously in 1988 by the 
Ginga LAC during a satellite manoeuvre (Makino et al. 
1988). The source had an average flux level of 26 mCrab 
(1-40 keV), and a power law spectrum with a — 1.7. Obser- 
vations both a month before and after by the Ginga ASM, 
by TTM in 1989 (In't Zand 1992) and by ROSAT in 1990 
and 1992 (Barret, Motch & Pietsch 1995) found comparable 
flux levels. Temporal analyses of both the Ginga detection 
and ROSAT data yielded a featureless / _1 power spectrum 
extending from 10~ 4 — 500 Hz (Tanaka & Lewin 1995; Bar- 
ret, Motch & Pietsch 1995), with neither QPO nor pulses 
being detected. 

Despite its detection by Ginga, the source had not been 
previously catalogued. Neither were X-ray bursts detected 
by Ginga. Together with its similarities to Cyg X-l (hard 
X-ray spectrum, strong flickering), this led to an early sug- 
gestion by Tanaka (1989) that it was a soft X-ray transient 
with a possible black-hole primary. Later, Strickman et al. 
(1996) called this suggestion into doubt, following examina- 
tion of data from CGRO/OSSE observations. They found 
that fitting both the Ginga and OSSE spectra produced a 
model with an exponentially cutoff power law plus reflection 
term. The observed cut-off energy of ~58 keV is typical of 
the cooler neutron star hard X-ray spectra. The recent re- 
port of three X-ray bursts detected by BeppoSAX (Ubertini 
et al. 1997) and our detection of optical bursts here confirms 
the presence of a neutron star accretor. 

Following the first ROSAT/PSPC all-sky survey ob- 



servations in September 1990, and the determination of a 
preliminary X-ray position, a search for the counterpart 
yielded a time variable, UV-excess, emission line star (Motch 
et al. 1994; Barret, Motch & Pietsch 1995). This source had 
B = 19.7 ± 0.1, and an uncertain V magnitude of V ~ 19.3, 
due to contamination by a nearby star. There was also evi- 
dence for ~ 0.3 m variations on a one hour timescale, but the 
time sampling was fairly poor. For this reason, we included 
this object in our target list for a high-speed photometry run 
at the South African Astronomical Observatory (SAAO). 
Further time-series photometry was also obtained on the 
William Herschel Telescope, La Palma (WHT) to confirm 
the variability that was seen. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 SAAO 1996 June 

Observations of a small (50x33 arcsecs) region surround- 
ing the optical counterpart to GS 1826-24 were made using 
the UCT-CCD fast photometer (O'Donoghue 1995), at the 
Cassegrain focus of the 1.9m telescope at SAAO, Suther- 
land on 1996 June 20. The UCT-CCD fast photometer is a 
Wright Camera 576x420 coated GEC CCD which was used 
here half-masked so as to operate in frame transfer mode, 
allowing exposures of as short as Is with no dead-time. We 
observed from UT 00:03:27 to UT 04:15:15, obtaining 1507 
consecutive 10s exposures. The conditions were generally 
good, being almost photometric, and the mean seeing was 
~1.5 arcsec. However, problems with telescope focus drift 



© 0000 RAS 



2 L. Homer et al. 



did lead to a 50 min interval of poorer quality data. In or- 
der to obtain the best signal-to-noise for the highest time 
resolution we used no filters. 

2.2 WHT 1997 July 

Observations of GS1826-24 were made on 1997 July 31 using 
the auxiliary port camera of the WHT. In order to provide 
the highest possible time resolution the TEK chip was win- 
dowed down to a small 33x33 arcsecs field-of-view, giving 
a dead-time of 10s. A series of 335 consecutive 25s B band 
exposures were taken from UT 21:54 to UT 01:10. Condi- 
tions were photometric throughout the run, and the seeing 
was generally quite good, varying between 1 and 1.7 arcsecs, 
apart from the first 15 min with ~ 2 arcsec. 

2.3 Photometry 

The imperfect auto-guiding of the SAAO 1.9m led to the 
positions of stars moving both randomly and with a definite 
trend during the run (from examination of coordinates of a 
bright star on the frame). To cope with the motion of the 
stars relative to each frame, an ARK routine (MAPCCD) 
was used to map the positions of about 10 bright stars on 
each image, compared to a reference image, and calculate 
initial offsets for star centres. A further improvement was 
then made by using these initial offsets to find the centroids 
of three bright stars, producing the final offsets. Similarly, 
two bright reference stars were used to calculate the small 
offsets relative to a reference image for the WHT run in one 
step. 

The following was then undertaken for both data sets: 

(i) The centres of all stars (above a chosen threshold) were 
measured on the reference image using the FIND routine in 
DAOPHOT II (Stetson 1987). The star list was edited to 
remove stars affected by any CCD defects as well as those 
stars both too close to the edge and very faint, creating a 
master star list. This included a suitable number of bright 
local standards as well as stars of comparable brightness to 
the counterpart. 

(ii) Due to moderate crowding of the counterpart with a 
nearby but fainter neighbour and the variable seeing, point 
spread function (PSF) fitting was essential to obtain good 
photometry. We performed the photometry reductions with 
the IRAF implementation of DAOPHOT II. An automated 
PSF fitting set-up was employed, whereby the same well- 
isolated and bright stars (see Fig. ^) were selected as PSF 
stars for each run, and the PSF was modelled on the profiles 
of these stars for each image separately. Then for each star 
of interest the PSF was scaled according to an initial aper- 
ture photometry estimate and fitted in order to measure its 
brightness. 

(iii) In order to double-check the psf fitting photome- 
try, we also ran aperture photometry using a rectangular 
aperture, positioned to include both the counterpart and 
its neighbour. The resulting lightcurve was identical, apart 
from being somewhat noisier. 

(iv) For this variability study only the relative brightness 
of a star is of importance, and so we applied differential 
photometry to help reduce the effects of any variations in 
transparency, employing up to 6 bright local standards. 



Owing to the variable seeing conditions and problems 
with telescope focus drift during the run, a small portion 
of the SAAO data is of poor quality. Hence, before subse- 
quent analysis those data points corresponding to frames 
with estimated seeing >2.5 arcsecs were deleted. The re- 
sulting 4.19 hour and 3.26 hour duration lightcurves for the 
optical counterpart to GS 1826-24 and a suitable star of com- 
parable brightness are shown in Fig.^, with a quadratic de- 
trend applied. Two features are immediately apparent; two 
very clear optical bursts, and evidence for short-term varia- 
tions in brightness. 



3 OPTICAL BURSTS 

A simple burst profile consisting of a linear rise, plus expo- 
nential decay was fitted to a 1200s section of each lightcurve 
including the bursts (as shown in Fig.^|). When fitting the 
burst a number of different models were used for the variable 
persistent flux; a double sinusoid constrained to fit the entire 
lightcurve and a linear or a quadratic baseline over the 1200s 
section, however the burst parameters were found to be in- 
sensitive to the model chosen. The time sampling (10 and 
25s at SAAO and WHT respectively) was adequate to con- 
strain the decay of the burst, but not the fast rise and hence 
peak flux. We can therefore only place limits. The fit to the 
SAAO June 1996 burst (reduced \ 2 = 1-6 for 113 d.o.f.) 
yielded the following parameters: rise time < 20s, e-folding 
decay time= 53.3 ± 2.7s, AF opt = HOOctss" 1 (i.e. 2.0 x av- 
erage flux). Similarly for the WHT August 1997 burst (re- 
duced x 2 = 0.8 for 29 d.o.f.), we obtained a rise time < 65s, 
e-folding decay time= 69.5±4.3s, AF opt = 390ctss _1 (note: 
constraining the rise time to < 20s gives AF opt = 500cts s _1 , 
i.e. 2.4 x average flux). The results are summarised in Table 

Many optical bursts have been detected from the X-ray 
burster X1636-536 (see e.g. Pedersen et al. 1982). In general, 
they exhibit a wide range of peak fluxes and profiles, with 
F m ax/ F pers = 1.4 — 4.8 ( in white light) and typical decay 
times ~few tens of seconds. Clearly, the bursts from GS1826- 
24 are not dissimilar. 

Lastly, we may compare the corresponding X-ray peak 
fluxes of these optical bursts to those of the bursts detected 
by BeppoSAX (Ubertini et al. 1997), using the results of 
Lawrence et al. 1983. They derived a simple power law rela- 
tion between the changes in the U, B and V band fluxes and 
corresponding X-ray flux variations during a well-studied 

burst of X1636-536, with = where 

a varies with the passband. 

The results of these calculations are presented in Table 
|i| where we use the a value for V as an approximation for 
the white light observations. The peak fluxes measured by 
BeppoSAX range from 350 to 520 mCrab, hence our lower 
limits of 190 and 270 mCrab imply consistency. 



4 SHORT PERIOD MODULATIONS 

Firstly, we detrended the extracted light curves for each of 
the stars using quadratic fits. In the case of GS 1826-24 the 
data points corresponding to the burst were also excluded. 
The maximum of the subtracted polynomial was constrained 
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Figure 1. Images showing the field of GS1826-24 (North up, East to the left): 10s white light exposure from SAAO June 1996 - 49x33 
arcsecs field of view (left), 25s B band exposure from WHT 1997 - 33x33 arcsecs field of view (right). The optical counterpart and the 
bright stars used for PSF fitting are labelled in each case. 
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Figure 2. Lightcurves of GS1826-24 and star of comparable brightness, a) High speed photometry in white light from SAAO with 10s 
resolution, b) conventional B band photometry from WHT with 35s time resolution. 



Table 1. Fitted parameters of a simple burst profile, with linear rise-time and exponential decay 



Date Time (UT) Rise time (s) AF op t(ctss 1 ) e-folding decay time (s) reduced x 2 

20 June 1996 02:18 < 20 1100 53.3 ±2.7 1.55 (113 d.o.f) 

1 August 1997 01:02 < 65 390 69.5 ±4.3 0.8 (29 d.o.f.) 

set < 20 500 69.1 ±2.8 0.9 (29 d.o.f.) 



© 0000 RAS, MNRAS 000, 000-000 



4 L. Homer et al. 



2500 

m 

m 
—> 

§ 2000 
o 

o 

la 

S 1500 
E 

■S iooo 



1000 



CO 

w 

4-J 

a 

o 
o 



800 



600 



400 



1 1 1 1 1 

a) SAA0 1996 June 20 
~ (White light) 


A 

\ 






/ 

| \ 




- 

A A A 

r 1 i t t i 


I A 
+ ) 

1 V 

1 c 

1 < 

1 ^ 

A A 
A 

1 1 


- 

- 

i i 


-500 





500 


Relative 


time (s) 




i i i i 

_ b) WHT 1997 August 1 
(r> banaj 

- 


i 

1, 

l\ 
1 \ 

i ^ 

' \ 


— ■ — ■ — i — n 
- 





1 \ 

t \ 

i \ 

i 

i V 





A_ A- A- A- *~ A- A- * ' — A- A- *- A ^ * 

1 1 1 1 1 1 


i 

i 


* ^ A * — A * ' ' ' ' . J 


-500 





500 


Relative 


time (s) 





Figure 3. Optical burst lightcurvcs (symbols) plus model fit (dashed line), a) High speed photometry in white light from SAAO with 
10s resolution, b) conventional B band photometry from WHT with 35s time resolution. The limiting model parameters are adopted 
with a rise time of 20s in each case. 



Table 2. Calculation of the peak X-ray burst flux corresponding to the two optical bursts of GS1826-24 



Date Passband a AF op t(ctss 1 ) ( ^° p '' l "°' j Fx.pe.rs (mCrab) Fx,maa;(mCrab) 



20 June 1996 White light 2.8 > 1100 2.0 27 > 190 

1 August 1997 B 2.55 > 500 2.35 31 > 270 



a Estimates of persistent flux taken from the XTE/ASM 2-10keV one-day average on the observation date 



to the time of least airmass for the field, ensuring reasonable 
compensation for the slow variations due to airmass changes. 
To search for periodic modulations two different methods 
were employed: (i) we calculated a Lomb-Scargle (LS) peri- 
odogram routine on each data set, to search for sinusoidal 
modulations (this periodogram is a modified discrete Fourier 
transform (DFT), with normalisations which are explicitly 
constructed for the general case of time sampling, including 
uneven sampling, see Scargle 1982); (ii) we constructed a 



phase dispersion minimisation periodogram (PDM), which 
works well even for highly non-sinusoidal light curves (see 
Stellingwerf 1978). 

The resultant LS periodograms for both GS1826-24 
datasets show distinct peaks at close to 11. 5d -1 (P = 
2.1 hr) (see Fig. ^). The exact measured frequencies are 
ll.0i0.2d" 1 and 12.0i0.4d" 1 for the SAAO and WHT 
lightcurves respectively, with semi-amplitudes of 0.05 m and 
0.07 m for a fitted sinusoid. Although the frequencies are only 
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Figure 4. Lomb-Scargle Pcriodograms (left) Phase-Dispersion Minimisation Pcriodograms (right) for GS1826-24 lightcurves from a) 
SAAO and b) WHT with bursts excluded. The PDM plots have been truncated at the lowest frequencies corresponding to the reciprocal 
of the data duration. 



marginally consistent (at ~ 2a level) within the formal ran- 
dom errors quoted, it must be noted that the datasets span 
only 2.0 and 1.6 cycles, with airmass changes of 1.0-2.4 and 
1.6-2.0 respectively, and hence the precision of the period 
determination is limited to ±0.8d _1 according to the exact 
form of detrending used. The formal confidence of these peak 
frequencies are greater than 99.98 per cent, as determined 
from a cumulative probability distribution (CDF) appropri- 
ate for the dataset (see Homer et at 1996 section 3 for details 
of the method employed). There are also significant peaks 
at the higher frequencies of 25.9 ± 2.3d" 1 , 33.4 ± O^d" 1 
(P = 55.7 ± 0.5 min, P = 43.1 ± 0.3 min) [SAAO] and 
21.4 ± O^.d" 1 (P = 67.3 ± 1.6 min) [WHT] correspond- 
ing to the shorter timescale variations. The results of PDMs 
confirm those of the Fourier analysis, with peaks occurring 
at 11. 5d~ x and ll.Od -1 , corresponding to the longer periods 
and also at 25.8d~\ 32.8d~ 1 and 21. 2d" 1 (see Fig. §). 

We also examined the LSP results for stars of compa- 
rable brightness. None of the four in the SAAO dataset ex- 
hibit any peak on 2-3 hr timescales with amplitude greater 
than 0.015" 1 , nor do these variations bear any phase relation- 
ship with the 2hr modulation seen in GS1826-24. Since the 
brightest stars in the field also show such small scale vari- 



ations, they are consistent with the small amplitude atmo- 
spheric transparency variations found at SAAO, enhanced 
by the effect of poorer photometry during the interval of 
inferior seeing. 

Similarly, the largest peak in the 2-3 hr range shown 
by the 3 check stars in the WHT dataset corresponds to a 
very low (0.007 m ) amplitude modulation. Again this is most 
probably atmospheric in origin. 

In any case, the presence of the 2hr modulation of 
GS 1826-24 in both datasets provides convincing evidence of 
its reality. Nevertheless, the short data span and the compli- 
cation of additional variability makes it impossible to con- 
strain the stability of the modulation and hence identify it 
categorically as orbital in origin. Longer time base observa- 
tions of this object are clearly necessary. 



5 DISCUSSION 

The period distribution of low-mass X-ray binaries (LMXB) 
has shown a scarcity of systems below 3 hr presumably in 
part due to their faintness, but there is a notable absence of 
any systems with periods between ~1 and 3 hr (White 1985; 
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Table 3. Properties of GS1826-24 and the short-period transient bursters" 



Source P a rb V B-V, U-B E B _ V Fx Active Quiescent 

(hr) ( M Jy) 

GS1826-24 2| 19.3 | 0.4 C , -0.5 C 0.4 C 30 c 1988-presentQ before 1988 

X0748-678 3.82 >23^-16.9 0.1, -0.9 0.42 0.1-60 1985-presentQ before 1985 

X2129+470 5.24 16.4-17.5 0.65,-0.3 0.5 9 <1979 - 1983-f 1983-present 

X1658-298 7.11 18.3 0.45,-0.4 0.3 <5-80 prc-1980s t 1980-present 



a Data taken from van Paradijs (1995) and references therein, unless cited elsewhere. 

b This work. 

c Motch et al. (1994). 

d Makino & theGinga team (1988). 

e Wade et al. (1985). 

f White, Nagase & Parmar (1995). 



White & Mason 1985). Recently, King & Kolb (1997) have 
investigated the formation of neutron star LMXBs. They 
found that in order to produce the relatively large fraction 
of soft X-ray transients in the < 1-2 day range, the secon- 
daries must have 1.3 M Q < M2 < 1.5 Mq at the onset of mass 
transfer and be significantly nuclear evolved (provided that 
the SN kick-velocity is small compared to the pre-SN orbital 
velocity). This mass range ensures that the mass transfer 
rates driven by angular momentum loss are below the criti- 
cal rate needed for SXT behaviour in the required fraction of 
neutron star LMXBs. The large initial masses then account 
for the rarity of any short P < 3 hr systems. If indeed the 2 
hr modulation of GS 1826-24 is confirmed as orbital in ori- 
gin, this would make it of great interest within the context 
of transient neutron star LMXB formation and evolution. 

However, in many respects GS1826-24 does show simi- 
larities to other neutron star binaries with comparable pe- 
riods. The X-ray bursters X1636-536 and X1735-444 ex- 
hibit optical modulations of a similar amplitude on their 
orbital periods (3.80 and 4.65 hr respectively), plus irregu- 
lar variability on somewhat shorter timescales (White, Na- 
gase & Parmar 1995 and references therein) , although unlike 
GS 1826-24 these systems are not transient on a timescale of 
decades. The optical emission of LMXBs is dominated by 
the reprocessed X-rays from the accretion disc and com- 
panion. The origin of the underlying sinusoidal modulation 
is interpreted as the varying contribution from the X-ray 
heated face of the companion star (vanParadijs & McClin- 
tock 1995). As for the irregular variability, this is probably 
caused by changes in the spatial distribution of the repro- 
cessing material in the accretion disc and/or fluctuations 
in the central X-ray luminosity, as suggested in the case 
of X1636-536 (vanParadijs et al. 1990). Moreover, GS1826- 
24 has a high Lx / 'L op t(~ 500), very similar to that of 
the compact 41 min binary X1627-673, but lower than the 
Lx/L pt ~ 700 of the 50 min binary X1916-053, which is 
a higher inclination dipping source (once again these are 
persistent sources). Since this ratio is related in part to 
the physical size of the system and the disc reprocessing 
area available, the similarity supports the hypothesis that 
GS 1826-24 is also a relatively compact system. 



Furthermore, the non-detection of GS 1826-24 prior 
to 1988 is characteristic of the observed variability of 
the transient bursting systems X0748-678, X1658-298 and 
X2129+470 (White, Nagase & Parmar 1995), which have 
either been detected over many years and then gone into 
quiescence for a similar period of time or vice-versa. These 
transients all have known periods in the range 3.82-7.1hrs 
(see Table |[). 

Clearly, further high-speed optical photometric moni- 
toring is required in order to confirm the stability of the 
2hr modulation and hence its orbital origin. With only a 
short .ROSvlTTightcurve published to date (Barret, Motch & 
Pietsch 1995) , confirmation might be possible from a longer 
X-ray observation. However, the low-amplitude of the ob- 
served modulation (0.06 m ) implies a low-inclination system 
(< 70°), and hence X-ray dipping behaviour is unlikely to 
be seen. 



6 ACKNOWLEDGMENTS 

We are grateful to Erik Kuulkers for comments on an earlier 
version of this paper. We thank Francois van Wyck (SAAO), 
Chris Benn and Peter Sorensen (La Palma) for their sup- 
port at the telescopes. L.H. acknowledges the support of a 
PPARC studentship. The WHT is operated on the island 
of La Palma by the Royal Greenwich Observatory in the 
Spanish Observatorio del Roque de Los Muchachos of the 
Institutio de Astrofi'sica dc Canarias. 



REFERENCES 

Barret D., Motch O, Pietsch W., 1995. Astr. Astrophys., 
303, 526. 

Homer L. et al., 1996. Mon. Not. R. astr. Soc, 282, L37. 
In't Zand J. J. M., 1992. PhD thesis, University of Utrecht. 
King A. R., Kolb U., 1997. Astrophys. J., 481, 918. 

© 0000 RAS, MNRAS 000, 000-000 



Lawrence A. et al., 1983. Astrophys. J., 271, 793. 

Makino F., the Ginga team, 1988. IAU Circ No 4653. 

Motch C. et al., 1994. IAU Circ No 6101. 

O'Donoghue D., 1995. Baltic Astronomy, 4, 519. 

Pedersen H. et al., 1982. Astrophys. J., 263, 340. 

Scargle J. D., 1982. Astrophys. J., 263, 835. 

Stellingwerf R. F., 1978. Astrophys. J., 224, 953. 

Stetson P. B., 1987. Pub. Astron. Soc. Pac, 99, 191. 

Strickman M., Skibo J., Purcell W., Barret D., Motch C, 
1996. Astr. Astrophys. Suppl. Ser., 120, 217. 

Tanaka Y., Lewin W. H. G., 1995. In: X-ray Binaries, 126, 
eds Lewin W. H. G., van Paradijs J., van den Heuvel 
E. P. J., Cambridge University Press. 

Tanaka Y., 1989. In: Two topics in X-ray astronomy, 23rd 
ESLAB symposium, 3, ESA:Paris, Bologna. 

Ubertini P. et al., 1997. IAU Circ. No 6611. 

van Paradijs J., 1995. In: X-ray Binaries, 536, eds Lewin 
W. H. G., van Paradijs J., van den Heuvel E. P. J., 
Cambridge University Press. 

van Paradijs J., McClintock J. E., 1995. In: X-ray Binaries, 
58, eds Lewin W. H. G., van Paradijs J., van den 
Heuvel E. P. J., Cambridge University Press. 

van Paradijs J. et al., 1990. Astr. Astrophys., 234, 181. 

Wade R. A., Quintana H., Home K., Marsh T. R., 1985. 
Pub. Astron. Soc. Pac, 97, 1092. 

White N. E., Mason K. O., 1985. Space Sci. Rev., 40, 167. 

White N. E., Nagase F., Parmar A. N., 1995. In: X-ray 
binaries, p. 1, eds Lewin W. H. G., van Paradijs J., 
van den Heuvel E. P. J., Cambridge University Press. 

White N. E., 1985. In: The Evolution of Galactic X-ray 
Binaries, 227, eds Triimper J., Lewin W. H. G., 
Brinkman W., NATO ASI, D. Reidel. 



Evidence for a 2 hr Modulation in GS1826-24 



© 0000 RAS, MNRAS 000, 000-000 



